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Abstract 

We propose a model of continuous opinion dynamics, where mutual interactions can be both positive and negative. 
Different types of distributions for the interactions, all characterized by a single parameter p denoting the fraction of 
negative interactions, are considered. Results from exact calculation of a discrete version and numerical simulations of 
the continuous version of the model indicate the existence of a universal continuous phase transition at p - pc below 
which a consensus is reached. Although the order-disorder transition is analogous to a ferromagnetic -paramagnetic 
phase transition with comparable critical exponents, the model is characterized by some distinctive features relevant 
to a social system. 
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1. Introduction 

Quantitative understanding of individual and social dynamics has been explored on a large scale 11 
in recent times. Social systems offer some of the richest complex dynamical systems, which can be studied using the 
standard tools of statistical physics. With the availability of data sets and records on the increase, microscopic models 
mimicking these systems help in understanding their underlying dynamics. On the other hand, some of these models 
exhibit novel critical behavior, enriching the theoretical aspect of these studies. 

Mathematical formulations of such social behavior have helped us to understand how global consensus (i.e.. 



agreement of opinions) emerges out of individual opinions [9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21.l 22ll23 



24 12511 . Opinions are usually modeled as variables, discrete or continuous, and are subject to spontaneous changes 
as well as changes due to binary interactions, global feedback and even external factors. Apart from the dynamics, 
the interest in these studies also lies in the distinct steady state properties: a phase characterized by individuals with 
widely different opinions and another phase with a major fraction of individuals with similar opinions. Often the 
phase transitions are driven by appropriate parameters of the model. 

In this paper we study a model of opinion dynamics by considering two-agent interactions. Continuous opinion 
dynamics has been studied for a long time ll26il27il28ll . with the models designed in such a way that eventually the 
opinions cluster around one (consensus), two (polarization) or many (fragmentation) values. The average opinion 
or macroscopic behavior have been emphasized only in some recent works [23, 24], where a phase transition from 
ordered to disordered phase has also been reported. However, in contrast to these models, we obtain here an ordered 
phase where even in the presence of a dominant opinion (symmetry broken phase), opposing opinions survive and 
a disordered phase where all opinion values coexist without any preference to any value (symmetric phase). Thus 
we present this in the general context of an order-disorder transition similar to that of the Ising and related models. 
We also compare our results with earlier works where a mean-field phase transition was observed in presence of 



contrarians in the society H21I1 . 

The paper is organized as follows: in Section 2 we introduce the model. Then in Section 3 the main results are 
presented along with the calculations and numerical simulations. In Section 4 we extend the model to include bond 
dilution and present the phase diagram. Finally we discuss our results in Section 5. 
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2. The Model 



We propose a new model for emergence of consensus. Let o,(f) be the opinion of an individual / at time t. In a 
system of individuals (referred to as the 'society' hereafter), opinions change out of pair-wise interactions: 

0,(f+l) = Oi(t)+ldijOj(t). (1) 

One considers a similar equation for Oj{t +1). The choice of pairs {/, j} is unrestricted, and hence our model is defined 
on a fully connected graph, or in other words, of infinite range. Note that this is simply a pair- wise interaction and 
we imply no sum over the index /'. Here fiij are real, and it is like an interaction parameter representing the influence 
of the individual with whom interaction is taking place. The opinions are bounded, i.e., -1 < o,(0 < 1. This bound, 
along with Eq. ([TJ defines the dynamics of the model. If, by following Eq. ([T]i the opinion value of an agent becomes 
higher (lower) than 4-1 (-1), then it is made equal to +1 (-1) to preserve this bound. The ordering in the system is 
measured by the quantity O - \ Oi\/N, the average opinion, which is the order parameter for the sy stem. 

The present model is similar in form to a class of simple models proposed recently (23l 13, Hs, 29 1, apparently 



inspired by the kinetic models of wealth exchange fs^j'sT]. A spontaneous symmetry breaking was observed in such 
models: in the symmetry broken phase, the average opinion is nonzero while in the symmetric phase, the opinions of 
all individuals are identically zero indicating a 'neutral state' . The parameters representing conviction (self interaction) 
and influence (mutual interaction) in these models were considered either uniform (a scalar) or in the generalized case 
different for each individual, i.e, given by the components of a vector. In addition to this there is an added feature of 
the randomness in the influence term which in eff'ect controls the sharpness of the phase transitions in these models. 

In our proposed model, the conviction parameter or self interaction parameter is set equal to unity so that in absence 
of interactions, opinions remain frozen. In such a situation, it has been observed previously that any interaction, 
however small, leads to a highly unrealistic state of all individuals having extreme identical opinions (either o,- = 1 V 
/ or Oj = - 1 V i) fl^ when the interactions take up positive values only. This suggests that one should generalize the 
interactions to include both positive and negative values. This is realistic also in the sense that it reflects the fact that 
in a social interaction of two individuals, there may be either agreement or disagreement of opinions. We therefore 
consider not only a distribution of the values of fijj (to maintain the stochastic nature of the interactions) but also allow 
fiij to have negative values. We define a parameter p as the fraction of values of pij which are negative, which, we will 
show later, leads to characteristic ordered and disordered states as in reality. 

The fact that we allow random positive and negative values for the interactions may suggest that the model is 
analogous to a dynamic spin glass model 1.32. ^3], as in the latter, one can consider a dynamic equation for the 
spins which formally resembles Eq. ([1}. However, the two dynamic models are not equivalent with the following 
differences: (i) the interactions in the opinion dynamics models are never considered simultaneously and thus the 
question of competition leading to the possibility of frustration does not arise, and (ii) there is also no energy function 
to minimize, (iii) the symmetry p — > 1 - p does not exist in our model, which is naturally present for spin-glass. We 
will get back to the comparison of the two models in the context of phase transition later in this paper 

The effect of negative interactions was considered previously in a different opinion dynamics model under the 
name Galam contrarian ll2l1l . The discrete, binary opinion model followed a deterministic evolution rule for a group 
of three or more individuals. It was shown that depending on the concentration of the contrarians, the system will 
either reach an ordered state, where there one of the opinions will have majority, or a disordered state, where no clear 
majority is observed. The critical behavior of the model is similar to the one we present here at least in the fully 
connected graph. However, our model considers continuous opinion values. Also, the Galam contrarians always take 
the opinion opposite to that of the majority. However, in our case we also consider the present state of opinion of the 
agents and accordingly even the discrete version of our model has three states. A two-state discrete version of this 
model will not show any ordered state. 



3. Results 

Unless otherwise mentioned, we keep pij values within the interval [-1,1] for simplicity. In principle, several 
forms can be considered for pij (annealed, quenched, symmetric, non-symmetric etc.). Further, there can be several 
distribution properties for pij in the interval [-1,1] (discrete, piecewise uniform and continuous distributions). Unless 
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otherwise stated, in our study, we would discuss the case when are annealed, i.e., they change with time. In other 
words, at each pairwise interaction, the value of fijj is randomly chosen respecting the fact that it is negative with 
probability p. For this case, the issue of symmetry does not arise. We consider distributions for both continuous and 
discrete yu,y. 

In all the above cases, we find a symmetry breaking transition. Below a particular value pc of the parameter p, the 
system orders (i.e., the order parameter O has a finite non-zero value), while the disordered phase (where (9=0) exists 
for higher values of p. Since this phase transition is very much like the thermally driven ferromagnetic-paramagnetic 
transition in magnetic systems, we have considered the scaling of the analogous static quantities, which are: 

(i) the average order parameter ((9), (. . .) denoting average over configurations, 

(ii) the fourth order Binder cumulant U - I - j^qt^, 

(iii) a quantity analogous to susceptibihty per spin, which we write as V = N[{0^) - (O)^]. 

We also calculate (iv) the condensate fraction fc - fi + / i, where /i and /_i denote the fraction of population 
with opinion 1 and -1 respectively, fc is exclusive to this class of opinion dynamics models and is expected to show 



Before a more detailed characterization of this transition, we find it useful to describe the behavior of the prob- 
ability distribution of the order parameter first. This distribution itself shows the signature of a phase transition. We 
consider both the cases with discrete (±1) and continuous (within the bound -1 < fi/j < +1) values for ju^. For 
each of the above case we consider both polarized and random initial conditions. We show the cases where the initial 
condition is fully polarized (all -i-l) and fijj are discrete (Fig. \Tla)) and continuous (Fig. [TJb)). We consider two p 
values, one below and one above the critical point. As can be seen from the figures, for p = 0.6 > pc the distributions 
are symmetric in both the cases, which manifestly imply disorder, while the ordered phase is asymmetric with a bias 
towards -i-l due to the initial condition. 

Now consider the case when the initial condition is random instead of being polarized as above. Then of course in 
a particular ordered state either majority positive or majority negative opinion can occur. Hence the order parameter 
distributions (after averages over many configurations) become symmetric (Figs[nc),(d)) even in the ordered phase 
for both discrete and continuous ju,y values. An interesting point, however, is the depletion of population in the 
intermediate (between the extremist with opinions +1) opinion values in an ordered society (FiglTJc)). This is even 
more pronounced when we make fi/j values discrete (±1). In this case, the distribution becomes discretized almost 
simultaneously with the symmetry breaking transition such that the intermediate opinion values do not exist at all in 
the ordered state (FigHfd)). Therefore, we find that in the ordered state with discrete fiij, the society becomes highly 
clustered in the sense that continuous variation of opinion is no longer a possibility. This of course is commensurate 
with the continuous fXij case described above, where possibility of ordering leaves the society highly clustered as is 
generally seen before and after decisive elections. 

It is important to note here that although the distribution of average opinion is sensitive to the initial condition, but 
the order parameter itself is not as the latter is obtained by taking the absolute value of the average opinion in a single 
configuration and then taking a further average over all configurations in a numerical study. Moreover, the critical 
behavior is also unaffected by the changes in initial condition. 

Our major interest lies in identifying the critical behavior in these models. Assuming that there exists a steady state 
for these models (which are numerically observed), one can derive exact expressions for the steady state probabilities 
/i, fo and /_i in the annealed discrete case where we assume the initial condition to be such that the agents have 
Oi(t = 0) G {- 1 , 0, H- 1 ). The scaling behavior of the order parameter and fc can be exactly obtained in this method. To 
do that, we consider the probabilities that an agent's opinion gets decreased (for initial opinion -i-l or 0) or increased 
(for initial opinion -1 or 0) or remains constant when two agents interact. For example, let us consider the change 
for the agent A, interacting with a second agent B (one can consider updated values of both but it does not matter): 
When both have opinion -i-l (probability of occurrence f^), A's opinion decreases with a probability p, giving the 
joint occurrence probability of the event as pf^. On the other hand, when B has opinion -1, the similar event has the 
probability (1 - p)f\f-\- One can similarly find out the probabilities of decrease and increase in all possible cases. 
The exact expression for the net increase probability is p f^^ + (1 - p)f\f) + pf)f-\ + (1 - P).f\.f-\ ™d that for the net 
decrease probability is pf^ + p fi fo -t- (1 - p)f(if-\ + (1 - p)f\f-\ - In the steady state these two should be equal, i.e.. 




Pfl + Pfifo + (1 - p)/o/-i = pf^i + (1 - p)fifo + pfof-u 



(2) 
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Figure 1: Probability distribution of opinions for polarized (all +1) initial condition with (a) discrete fijj {pc = 1/4) and (b) continuous t^jj 
(Pc ~ 0.34); same for random (uniform between [~1 : +1]) initial condition with (c) discrete fiij {pc = 1/4) and (d) continuous fijj {pc « 0.34). All 
data are for N = 4096. 



which simplifies to 

(2/i+/„-l)[p-/o(l-p)]=0. (3) 
This means, either 2/i + /o = 1, i.e., /i = (1 - /o)/2 = f-i which implies a disordered phase, or 

/o = T-^. (4) 

1 -p 

Next we show that at criticality all three fractions would become equal such that Pc = I /4. Let us take the solution in 
the disordered phase where fi - /_i . We consider the processes contributing to the in/out flux for /q. We enumerate 
all possibilities as before and get the following: flux into /o is 2[(1 - /o)/2]^ and flux out of fo is /o(l - /o). So, in the 
steady state 

''^f = ^. (5) 

Hence, either /o = 1, which can be ignored by considering the steady states of the other two fractions, or /o = 1 /3. 
Now, just at the critical point this solution will begin to be vaUd. Hence, at critical point all three fractions are 1 /3, 
leading to = 1/4. 

In the ordered phase, the order parameter O is given by |/i - /- 1 1 . To evaluate /i and /_ i , we calculate the flux in 
and out of /i . Flux out of /i is pff + (1 - [1 - (/i + fo)] and flux into /i is (1 - /?)/o/i + /?/o[l - (/i + fo)]- Hence 
at steady state, /i is given by 

^' " 2(l-2p + p^) ' 

where we have used Eq. (|4]l. Hence the order parameter is given by 



^ ^ l-3p + 2p^± Vl-6p + V_4p3 ^ 2p-l 

(1-/7)2 l_p- 

The variation of the order parameter with p is shown in Fig. [2] It shows the expected behavior, i.e., it vanishes at 
p = Pc- Rewriting the above equation in terms of x = p^ - p, algebraic simplifications give 

_ 3/8 + 2x + 2x2 + ^9^/4 _ 3^2 2x+l/2 
^ ~ 9/16 + 3.=c/2 + x2 3/4 + x ■ 

As X — > 0, (9 ~ y/x, implying the critical exponent for O is /3 - 1/2. This also agrees well with the power law fit of 
the order parameter expression (Eq. (|8)) near the critical point (Fig. Calculation for f on the other hand shows 
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Pc-P 

Figure 2: Discrete power law beliavior of the order parameter O near tlie critical point pc (Eq. (8)) showing /? = 1/2. The dotted line is ;c''-. 
a visual guide. Inset: (a) Phase diagram. The points represent simulation results. They are in better agreement with the analytical results as the 
system size is increased from 256 to 2048. (The lower half of the phase diagram follows from symmetry.) (b) Linear scaling of ~ /?■ The dotted 
line is x' . 

that it has a constant value fj^ -Ijli beyond pc and for p < p^, fc - fi^ ~ Jc', i.e., vanishes Hnearly at pc, which also 
perfectly agrees with numerical simulations. Numerical simulations for the model with continuous fiij yields similar 
behavior with fj^ ^ 0.22 since opinion values other than ±1 and exist. 

Using the same kind of argument as above, one can show that there will be no phase transition when yU,j = ±/vo 
with > 2. It is obvious that in this case the opinions can have only two values +1 and -1. Let /i be the fraction of 
agents having opinion +1. Once again considering net decreases and increases of opinions and equating them in the 
steady state, we get 



(9) 



For any non-zero p the solution of this equation is /i - 1/2 thereby giving complete disorder Hence, in this condition 
there cannot be an ordered phase for any finite p. 

We performed Monte Carlo simulation for different system sizes (A^ = 64,256,512, 1024,2048) to estimate pc 
and all the relevant exponents for discrete as well as continuous ////s (see Fig. [3] where the data for the continuous 
distribution of yu,y are presented). A Monte Carlo step is the simultaneous update of agent's opinion values. For each 
simulation point, sufficient relaxation time was given (depending on system size), such that the measurable quantity 
reached a steady-state value. Then the ensemble average of those steady state values were taken (number of ensemble 
again depends on system size, see figure captions for details). 

We estimated the critical point pc from the crossing of the Binder cumulants for different system sizes 13411 . Our 
estimate is pc - 0.249 + 0.001 for the discrete case which is consistent with the analytical value of 1/4 derived 
earlier The critical Binder cumulant is U* = 0.30 ± 0.01. For the continuous case, pc - 0.3404 + 0.0002, and 
U* = 0.284 + 0.004. We find excellent finite size collapse for all cases. We estimate the correlation length exponent 
V = 2.00 ± 0.01, the order parameter exponent p = 0.50 ± 0.01 and the fluctuation exponent y = 1.00 + 0.05. 

One can consider several other distributions with a similar parameter p to test the universality of the phase transi- 
tion: e.g., we take yU/y - 1 with probability (1 - p) and = -1/2 with probability p. The transition point shifts but 
the exponents remain same. 

We also conducted detailed simulations for the case where /i^ are quenched variables, i.e., they do not change 
with time. For the fully connected graph considered here, numerical results for continuous /i,/s (see Fig.|4]i indicate 
that critical behavior is the same as that of the annealed case mentioned above. 



4. Model with bond dilution 

Here we briefly consider the case when all the interactions between the agents are not allowed, i.e., agents are 
selective in their interactions. In the case of continuous distribution of it does not make sense to parametrize the 
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Figure 3: Data for continuous, annealed /j; j model, showing (a) finite size scaling of the Binder cumulant U for dilferent system sizes N; the critical 
point is pc = 0.3404 ± 0.0002, and the best data collapse is for v = 2.00 ± 0.01. Inset: Variation of U with p for different system sizes; (b) finite 
size scaling of order parameter O for different N; best data collapse is for /} = 0.50 ± 0.01. Inset: Variation of the order parameter O with p for 
different system sizes; (c) finite size scaling of V for dift'erent N; best data collapse is for y = 1.00 ± 0.05. Inset: Variation of V with p for dift'erent 
N. Number of averages for different system sizes are 3000 for N = 256, 1800 for W = 512, 1000 for N = 1024 and 400 for N = 2048. 



probability of jt/,j = 0. It is therefore useful to consider the discrete distribution only for the dilute case by considering 
fijj - with probability q along with yU/y = -1 with probability p as before and fitj - 1 otherwise. Thus, q is naturally 
the 'bond dilution' parameter 

Considering the balance of the increases and decreases in the order parameter in the steady state, one gets 

fiP + fofi p+fof-i[l-(p + q)] 

= A+fo.n[i-(p + q)]+fof-^p. (10) 

This gives, either /i = i.e., disorder, or in the ordered state 

/o = ^-^ . (11) 

Now, in and out fluxes of /o gives 

fip + 2/i/_i[l-(/, + ^)]+/_V=/o/i[l-(p + ?)] 

+ foflP + fof-l[l-(p + q)]+fof-lp. (12) 

In the disordered state, the only feasible solution is /o = |. As before, assuming the continuity of /o across the 
transition point, we get the phase boundary equation as 

Pc^\(l-qc), (13) 
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Figure 4: Data for continuous, quenched fjjj model, sliowing (a) finite size scaling of the Binder cumulant U for different system sizes N; the 
critical point is = 0.34 ± 0.01, and the best data collapse is for v = 2.00 ± 0.01. Inset: Variation of U with p for different system sizes; (b) finite 
size scaling of order parameter O for different N; best data collapse is for /? = 0.50 ± 0.01. Inset: Variation of the order parameter O with p for 
different system sizes; (c) finite size scaling of V for dift'erent N; best data collapse is for y = 1.00 ± 0.05. Inset: Variation of V with p for dift'erent 
W.Number of averages for different system sizes are 3000 for N = 256, 1800 for N = 512, 1000 for N = 1024 and 400 for N = 2048. 



where qc is the critical value of q on the phase boundary. This phase boundary is shown in Fig. |5] Clearly ^ = 
limit is the unrestricted case studied in the previous section. The fact that the ordered phase extends to = 1 is for 
p - Ois consistent with the fact that for an infinite dimensional lattice, percolation transition takes place atq = 1. The 
universality of the model along the phase boundary needs to be done. One can expect it to be same as before. 



5. Discussions 

We present a simplified model for opinion formation in a society of highly connected individuals. With the 
dynamical rule defined by Eq. ([TJ, an agent modifies his/her opinion under the influence of another This evolution 
rule corresponds to cases where opinions are modified following discussions/debates with another individual. The key 
feature of our model is the inclusion of negative interactions with probability p. We study the steady-state collective 
behavior of this model. While modeling the opinion formation, the effects of 'self-conviction' and 'external-pressure' 
have been considered before (see e.g., (|35[). 

Here we concentrate upon the study of 'consensus formation' as a dynamical critical phenomenon. Eq. ([T]i of 
course does not include all social complexities involved in such interactions, however in this simple form it manifests 
some intriguing features. 

With the introduction of a single parameter p, defined in a simple manner, the proposed model shows the existence 
of a universal phase transition and also some additional desirable features representing a real society. The parameter 
p plays a role of the disordering field (similar to temperature in thermally driven phase transitions). Beyond a certain 
value of the fraction of negative interaction, a phase transition from an 'ordered' state (with most of the individuals 
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Figure 5: Pha.se diagram for tlie dilute model. 



having opinions of same sign) to a 'disordered' state (where the opinions can have different signs and add up to zero) 
occurs (see Eq. (O). 

The highHght of our model (with discrete values of yU,j) lies in the unique selection of an ordered state with 
discrete opinion values (±1 and 0) while in the disordered state, all opinion values in [-1, +1] coexist (Fig. [He)). The 
disordered state is one with a lot of disagreement, hence all types of opinions co-exist in the society. But as it starts to 
get ordered (below pc), polarization occurs, and marginal opinions cease to exist, resembling the ordering in a multi- 
party election scenario. This unique selection of the ordered state is in fact independent of the initial conditions which 
adds to the richness of the model. This of course also happens in the bounded confidence models where opinions, 
originally varying continuously, ultimately cluster around a few values typically. However, as akeady mentioned, 
in these models, the dynamics is designed to achieve this, whereas in our model, it happens naturally without any 
imposed restrictions. 

The phase transition in our model presents a case of a classical continuous phase transition with simple exponent 
values and showing finite size scaling behavior To model 'social temperature' which essentially destroys consensus 
rather than forming it, we have kept negative interactions among the agents. This leads to the desirable feature that 
even in the 'ordered' phase, opinions of both signs can coexist and that the disordered state is also not a 'neutral state' . 
In contrast, in the models with only positive interactions 11231124 12511 . where each 'scattering' leaves the two agents 
closer to each others' opinion, finite size behavior was absent and the order of phase transition difficult to identify. 
Thus, we believe our model is closer to reality and also in terms of critical behavior presents significant modifications. 

Another important point to be mentioned is that one could change the distribution of the interactions in a way 
such that no 'neutral' individual remains in the society (opinion values are +1). It is intriguing and to some extent 
counter intuitive that such processes lead to no consensus in our model (see Eq. (jSj). When opinions can take values 
equal to ±1 only, a comparison with a Ising spin model is bound to arise. The model with which one should compare 
is the fully connected Ising model with random bimodal distribution of the interactions (-/ with probability p and 
-1-7 with probability (1 - p)). However, in that model there is indeed a phase transition occurring at p - 1/2 llssll 
from a ferromagnetic to a spin glass phase. On the other hand, we get a transition only when the opinions can take 
more values than just +1 (e.g., ±1 and 0). In our opinion dynamics model, the ordered phase may be regarded as a 
ferromagnetic phase and the disordered phase as a paramagnetic phase (certainly not a spin glass phase as the opinions 
do not attain a frozen state had it been so). 

When agents are 'selective' in their interactions and some of the interactions are absent or muted, this effectively 
makes those p/j = 0. If pij - with probability q, then we can conceive a phase diagram with respect to the parameters 
p and q (see Fig. [Sj separating the ordered and disordered phases. 

A comparison of our model with the model of Ref |21 ] which introduces Galam contrarians (see also ll36ll37[ 13811 ) 
may also be made. In the latter, one has two discrete choices of opinion while the general case of the present model 
involves continuous opinion values. The mean-field critical behavior, however, is observed in both models with same 
exponent values. In Ref 12 ill the contrarians would take the opinion which is exactly opposite of that of the majority 
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of the group. The evolution of opinions is clearly different in our case, where the original opinion of an agent is also 
considered while assigning the changed value. This makes even the discrete version of the present model consist of 
three states (-1, +1). Furthermore, we have shown that (see Eq. (|9}) if we consider the discrete version to have 
two states only, then there is no ordered state in the model as soon as any finite fraction of negative interactions is 
introduced. 

Results obtained from opinion dynamics models can be compared to real data to a certain extent. Of course, the 
microscopic rules governing the dynamics in a model cannot be verified directly but one can justify the model by 
comparing its macroscopic behavior with real data. Discrete opinion models may be applicable to election scenarios 
ll21l1 . Continuous opinion models, on the other hand, mimic the case of rating or degree of support for an issue. 
So in our model, -1 represents extreme unfavorable opinion, +1 is extreme favorable while zero means an average 
rating/indifferent response. Thus the order parameter in the model corresponds to the overall rating and ordered 
state means there is a clear-cut decision made. A disordered state means the absence of a decision. The situation is 
comparable to a case of final verdict arrived at by a panel of juries. If there is a lot of disagreement among the juries a 
decision is hard to achieve - this result is indeed obtained in our model where the disagreement is represented by the 
negative interaction terms. 

Although we have considered continuous opinions in the model, it was shown exactly that using discrete opinions 
1,-1,0 also leads to the same critical behavior From the statistical physics viewpoint, the important role in the phase 
transition is thus played by the parameter p quantifying the fraction of negative interactions, irrespective of the fact 
that opinions are continuous or discrete. This precisely indicates that p is the relevant parameter in the model and the 
nature of opinion is irrelevant as far as critical behavior is concerned 

We conclude with the remark that the values of the exponents yS, y are very similar to that of the mean field 
exponents of the Ising model. Interpreting v as v'd where d is the effective dimension in this long ranged model and 
putting c/ = 4 as is done in small world like networks (3^, the value of the effective correlation length exponent 
v' = 1 /2 also coincides with the mean field value. 

For future study, the properties of this model in various lattices and networks and also its dynamical behavior 
would be interesting [40]. 
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